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In a wide variety of materials, such as layered 
copper oxides, heavy fermions, organic salts, and 
the recently discovered iron pnictides, supercon- 
ductivity is found in close proximity to a mag- 
netically ordered state^ ^. The character of the 
proximate magnetic phase is thus believed to be 
crucial for understanding the differences between 
the various families of unconventional supercon- 
ductors and the mechanism of superconductivity. 
Unlike the antiferromagnetic order in cuprates, 
which is well described by the spin Heisenberg 
model, the nature of the magnetism and of the 
underlying electronic state in the iron pnictide 
superconductors is not well understood. Neither 
density functional theory nor models based on 
atomic physics and superexchange, account for 
the small size of the magnetic moment^. Many 
low energy probes such as transport '^, scanning 
tunneling microscop^^ and angle-resolved pho- 
toemission spectroscopy (ARPES)[6j measured 
strong anisotropy of the electronic states akin 
to the nematic order in a liquid crystal, but 
there is no consensus on its physical origin, and 
a three dimensional picture of electronic states 
and its relations to the optical conductivity in 
the magnetic state is lacking. Using a first prin- 
ciples approach, we obtained the experimentally 
observed magnetic moment, optical conductiv- 
ity, and the anisotropy of the electronic states. 
The theory connects ARPES, which measures one 
particle electronic states, optical spectroscopy, 
probing the particle hole excitations of the solid 
and neutron scattering which measures the mag- 
netic moment. We predict a manifestation of 
the anisotropy in the optical conductivity, and 
we show that the magnetic phase arises from the 
paramagnetic phase by a large gain of the Hund's 
rule coupling energy and a smaller loss of kinetic 
energy, indicating that iron pnictides represent a 
new class of compounds where the nature of mag- 
netism is intermediate between the spin density 
wave of almost independent particles, and the an- 
tiferromagnetic state of local moments. 

Below the Neel temperature of the order of 150 K the 
parent compounds of the iron pnictide superconductors 
remain metahic with a magnetization density oscihating 
in space (spin density wave, SDW). The sublattice mag- 
netization is concentrated on iron atoms and its arrange- 
ment in space is antiferromagnetic in the x direction and 
ferromagnetic in the y direction^. 

We use the combination of the density functional the- 



ory (DFT) with local density approximation (LDA) and 
the dynamical mean field theory (DMFT) |8] to study the 
archetypical iron pnictide compound BaFe2As2 in both 
the magnetic SDW and the paramagnetic (PM) state. 
The size of the theoretical magnetic moment is 0.86 /i^, 
similar to the measured moment of 0.87 /iBp^, but much 
smaller than 1.75 fisi^^ obtained by the local spin den- 
sity approximation (LSDA) within DFT. The consider- 
ably smaller magnetic moment obtained here is due to the 
fact that the competing PM metallic state is a correlated 
metal, which contains very fast fiuctuating moments in 
time, but no static moment. Only a small part of these 
fiuctuating moments acquires a static component in the 
ordered state. 

The onset of magnetic order has a profound impact on 
the electronic structure, and these changes are probed 
by optical spectroscopy. Figure flla) shows the in-plane 
(averaged over x and y direction) optical conductivity of 
BaFe2As2 in the SDW and PM states calculated by both 
LDA+DMFT and L(S)DA. Fig. [ijb) reproduces mea- 
sured in-plane optical conductivity from Refs. 10 andfTTl 
Both theory and experiment j^^ l "*^^ * show a reduction of the 
low frequency Drude peak, which indicates a removal of 
a large fraction of carriers in the ordered state. Our cal- 
culation captures all the important qualitative features 
measured in experiments. In both PM and SDW states 
there is a broad peak due to interband transitions cen- 
tered around 5500cm~^[12 . Below 2000 cm~^ the op- 
tical conductivity of the SDW phase shows a few extra 
excitations appearing as peak and shoulders centered at 
1250 cm~^ (arrow 2 in cyan), shoulder structure at about 
800 cm~^ (arrow 1 in blue), and a small peak at 1800 
cm~^ (arrow 3 in green). These additional excitations 
appear in experiment at slightly smaller energies, as seen 
in Figs.[TJa) and (b). 

These extra peaks and shoulders strongly depend on 
the polarization of the light, as shown in Fig. flic), where 
we plot separately the x and the y component of the 
optical conductivity. The first two excitations (blue ar- 
row 1 and cyan arrow 2) are much more pronounced in 
the X direction, while the third peak (green arrow 3) is 
more pronounced in the y direction. Also the conductiv- 
ity is considerably larger in the x direction (antiparallel 
spins) than in the y direction (parallel spins). At low 
frequency, the optical conductivity of the SDW phase 
shows Drude like behavior, with Drude weight consider- 
able smaller than in the PM phase. The theoretical value 
of the plasma frequency for the x, y and z direction (1.14, 
0.88 and 0.69 eV, respectively) are much smaller than 
the plasma frequency of 1.60 eV of the PM DMFT calcu- 
lation^^, which agrees very well with experimental esti- 
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FIG. 1: Optical conductivity and DOS of BaFe2As2. 
(a) calculated in-plane average optical conductivity by 
LDA+DMFT and L(S)DA in the SDW and PM states; 
(b) experimental in-plane optical conductivity in the SDW and 
PM states. Data are taken from Refs. 10 and [XT'; (c)the xx, 
yy, and in-plane average optical conductivity in the SDW 
state calculated by LDA+DMFT and LSDA (only the aver- 
age is shown for LSDA); (d)total DOS in the SDW and PM 
states calculated by LDA+DMFT and L(S)DA and (e)the 
projected DOS of Fe 3d xz, yz and xy orbitals in the SDW 
state calculated by LDA+DMFT, plotting positive/negative 
for majority/minority spin. 



mates^^ESl. This reduction was also observed experimen- 
tally^. The anisotropy of the in-plane optical conductiv- 
ity has not been studied experimentally, because of the 
fairly coarse spatial resolution imposed by the diffraction 
limit, spanning multiple magnetic domains. Near field 
optics^^, or techniques to prepare monodomain samples^, 
are promising avenues to test our prediction. 

It is useful to analyze the optical conductivity at var- 
ious frequency scales. The Drude weight is controlled 
by the Fermi surface size, and by the mass enhance- 
ment of the low energy quasiparticles. In the SDW 
state, the mass enhancement is smaller (2, 1.7, 1.7 and 
1.5 for the t2^/majority, t2^/minority, e^/majority, and 



e^f /minority orbitals) than in the PM state ( 3 and 2, for 
the t2g and eg orbitals, respectively). While the quasi- 
particles become lighter in the SDW phase, the Fermi 
surface area is much smaller, and the latter effect domi- 
nates, resulting in a reduction of the Drude weight. 

We also integrated the in-plane optical conductivity to 
obtain the effective kinetic energy of a low energy model, 
in both the SDW phase and the PM phase. At very 
low energies, the onset of magnetism results in an in- 
crease of the optical conductivity, due to the coherence- 
incoherence crossover. The long range order makes the 
material more coherent. Consequently the very low en- 
ergy model gains kinetic energy in the ordered state. At 
intermediate energies, however, kinetic energy is lost as 
the result of the opening of the SDW gap on the Fermi 
surface. 

Since the optical conductivity of the PM phase is quite 
temperature dependent, we compare the SDW phase 
and PM phase at the same temperature (T = 72. 5K), 
the latter being a metastable state at low temperature. 
Using this procedure, we find that the missing weight 
from opening the SDW gap is recovered only around 
10000 cm~^, many times larger scale than the gap value. 
This should be contrasted with the classic weakly cor- 
related materials, where the spectral weight is recovered 
immediately above the SDW gap. 

Recent dc conductivity and scanning tunneling mi- 
croscopy (STM) measurements have detected a large 
anisotropy in the ab plane.^^. The theoretical zero fre- 
quency limit of the optical conductivity is 7350, 4400, 
and 3350 Q~^cm~^ for the xx^ yy and zz components, 
respectively. This implies that the resistivity in the y di- 
rection is about 1.67 times of that in the x direction, in 
good agreement with the results of resistivity measure- 
ments of Chu et al.^. 

We now turn to electronic density of states. FigurJlld) 
shows the total density of states (DOS) in both the PM 
and SDW phases. In the SDW phase, DMFT DOS 
shows a clear pseudogap on the scale of 0.15 eV around 
the Fermi level, in good agreement with STM measure- 
ments^. The LSDA DOS also shows a pseudogap at the 
Fermi leveP^, however, its width is more than 0.5 eV, 
therefore LSDA misses the structure below 4000 cm~^ in 
the optical conductivity. 

The onset of stripe magnetic phase is also accompa- 
nied by a rearrangement of the iron crystal field states, 
which gives rise to orbital polarization. This polariza- 
tion is uniform in space (ferro-orbital ordering), as sur- 
mised by Singh^^. The partial density of states of an 
Fe atom is shown in Fig {lie). The minority density is 
given a negative sign. To extract the anisotropy of the 
electronic structure, we integrate the partial density of 
states ofxz and yz orbital {Axz{(^) and Ayz{uj)) to obtain 
their occupation, and evaluate their difference An (A) = 
J^^[A,,{uj) - Ay,{uj)]duj/{lJ^^[A,,{uj) + Ay,{uj)]duj). 
This defines the energy dependent orbital polarization. 
For large cutoff A, the orbital polarization is close to 
zero for majority electrons and around 0.13 for minority 



electrons. At low energy, the anisotropy is enhanced to 
1.23 (0.45) for majority (minority) carriers, when A is 
0.15 eV, the size of the optical SDW gap. 

The anisotropy of the partial density of states provides 
a natural explanation for the anisotropy in the optical 
conductivity; the yz density of states has less electronic 
states at the Fermi level and the main peak of the yz 
orbital is further away from the Fermi level compared to 
xz orbital. The optical conductivity in x direction comes 
primarily from xz and xy orbit als, and is thus larger than 
the conductivity in y direction, which is connected to yz 
and xy orbit als. 

In FigJ2|a), we show LDA+DMFT momentum- 
resolved electronic spectra A(k, cj) in the SDW phase. 
This quantity is probed by angle-resolved photoemission 
spectroscopy (ARPES). FigsJ2Vb) and (c) compare the 
Fermi surface of the PM and SDW states, displayed in 
the PM Brillouin zone. 
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FIG. 2: ARPES and Fermi surface of BaFe2As2. (a) A(k,a;) 
in the F plane in the path M' ^ F ^ M ^ F of the SDW 
state. (See the locations of M and M' points in panel (b)). In 
the path M ^ F we overlayed ARPES data from Ref.[6 . (b) 
Fermi surface in the PM state, (c) Fermi surface in the SDW 
state, plotted in the PM Brillouin zone. Blue arrows mark the 
position of the Dirac cones. Note that the reciprocal vectors 
in the order state are (1/2,-1/2,0) in the direction of the 
ferromagnetic ordering, and (1/2, 1/2, 1/2) in the direction of 
the antiferromagnetic ordering, hence points F and M' are 
equivalent, while points F and M are not. (d) A{k^uj) in the 
SDW state with shadow bands plotted by equal intensity for 
clarity. Arrows mark the three types of optical transitions 
which give rise to the three peaks in the optical conductivity. 

In the PM state, the topology of the Fermi surface is 
very similar to LDA predictions^^ with three cylinders 
centered at F point and two at M point. In the SDW 
phase, the Fermi surface of the LSDA calculation (not 



shown) and LDA+DMFT is very different. The mag- 
netic order reconstructs the Fermi surface into smaller 
more three dimensional pockets. Out of three cylinders 
centered at F point, one remains in the SDW phase. This 
cylinder does not intersect the F plane within LSDA (not 
shown), but has been clearly identified in experiment. 
The other cylinders reconstruct into more three dimen- 
sional pockets. 

In the SDW state, there are two inequivalent directions 
between F and M, here named M and M', pointing along 
the antiferromagnetic and ferromagnetic direction of the 
Fe-Fe bonds, respectively, (see Fig.|2Vb)). Graphene-like 
Dirac points were recently identified by ARPES along 
the antiferromagnetic directioiP. FigJ2[a) shows that a 
crossing of two bands occurs very near the Fermi level 
between F and M, at 3/4 of the way, marked by a white 
arrow. The crossing is below the Fermi level, hence the 
pocket is electron like. In Fig. |2|c) we mark the same 
tiny pocket by blue arrows, to show that it has indeed 
a shape of a Dirac cone. There are two such symmetry 
related Dirac cones in the F plane and two in the Z plane. 
Notice that these cones appear only in the path between 
F to M (antiferromagnetic direction) and not in the F to 
M'. 

In the most right-hand part of Fig.|2|a), we overlay our 
results with ARPES measurements of Richard et a/.^, to 
emphasize common features. The overall position of the 
bands is in very good agreement without any need of shift 
of the Fermi level or renormalization of the bandwidth, in 
contrast to common need for shifts and renormalization 
when comparing DFT-derived bands with ARPES. 

The LDA+DMFT Fermi surface also has good agree- 
ment with the ARPES measurement in the Z plane by 
Shimojima et al}^ . In particular, the red electron pock- 
ets centered at Z, which have a two fold symmetry and 
mostly xz character, were identified in Ref.18, 

In Fig. [2|d) we replot the momentum resolved elec- 
tronic spectra A(k, cj) without the SDW coherence fac- 
tors, to enhance the shadow bands. The arrows in this 
graph connect the features in the electronic structure 
with the peaks in the optical conductivity. We mark tree 
types of vertical transitions corresponding to the three 
peaks in Fig. [l] The first shoulder comes primarily from 
transitions within the xz orbital, namely between the fiat 
band around — O.leV and the hole pocket at F and M. 
These transitions are between a shadow band and a non- 
shadow band, hence they appear only in the SDW phase. 
The second peak in optics, marked by cyan arrow, comes 
primarily from transitions between the xz and xy or- 
bitals, with some transitions between non-shadow bands 
only, visible also in the PM state, and some additional 
transitions between a shadow and non-shadow band. Fi- 
nally, the third peak comes mostly from transitions be- 
tween the xy and yz orbitals, and mostly from transitions 
between a shadow to non-shadow band. 

In correlated materials new physics, such as supercon- 
ductivity, spin and orbital polarization, emerge from the 
competition between Coulomb interaction and kinetic en- 



ergy. 

A unique physical characteristics of iron arsenic mate- 
rials is that the kinetic energy loss in the SDW phase is 
compensated by a gain in Hund's rule coupling energy. 
Comparison of the SDW and PM histograms, describ- 
ing the probability of different iron configuration in the 
solid, shows that in the SDW state the high spin states 
become more probable. This results in an overall gain 
of the Hund's rule coupling energy of about 500 K per 
Fe. Overall kinetic energy is lost in the SDW state for 
about 300 K per Fe, resulting in a net energy gain of 
about 200 K per Fe. This is different from classical SDW 
transition where kinetic energy is compensated by the 
reduction of the Hubbard correlations. 

The competition between kinetic and correlation en- 
ergy takes different forms at different energy scales, and 
results in an unusual energy dependence of the spin 
and orbital polarization. Spin polarization affects most 
strongly the electrons far below the Fermi level. For ex- 
ample, the exchange splitting - as determined from the 
frequency dependent potential (self-energy) at high fre- 
quency - is three times larger than at zero frequency. This 
high frequency regime, is governed by the strong Hund's 
rule coupling on iron atom, enhancing magnetic moment. 

At low energies, in the SDW state, well defined quasi- 
particles form, and the residual Hund's rule coupling be- 
tween these quasiparticles is weak. To minimize the ki- 
netic energy loss in the SDW phase, the quasiparticles 
propagate mainly along the antiferromagnetic x direc- 
tion, the direction which is not blocked by the Pauli ex- 
clusion principle. This generates strong orbital polar- 
ization, but only at low energy, where the quasiparticles 
are well formed, and the effective Hund's rule coupling is 
weakest. On the contrary, the overall orbital polarization 
is weak. 

Our finding that spin polarization is larger at high en- 
ergy while the orbital polarization is most pronounced 



at low energies leads to definite predictions for the 
anisotropy of the optical conductivity and can be tested 
also by STM. 

Energy dependent polarizations and the enhancement 
of coherence of the low energy quasiparticles in the SDW 
phase, can only be described by frequency dependent po- 
tentials and Weiss fields, as is done in DMFT. This ex- 
plains the failure of static mean field theories such as 
LDA to capture both the correct moment, which lives at 
high energies, and the low energy spectra. In a renor- 
malization group picture of this material, one observes a 
different strength of the Hund's rule coupling at different 
energy scales. At high energy, Hund's rule coupling is 
very strong, while it fades away at low energy, but gives 
an imprint on the massive and anisotropic low energy 
quasiparticles. This is central for a proper description 
of the magnetic phase, and is likely to be important for 
the mechanism of the unconventional superconductivity 
in these materials. 

METHOD 

To show that the origin of the anisotropy is electronic 
rather than structural, we use the experimental lat- 
tice constants and internal coordinates of the param- 
agnetic phase^. We use the continuous time quantum 
Monte Carlo as the impurity solver and charge self- 
consistent version of LDA+DMFT, described in detail 
in ReflTOl We use ab initio determined Coulomb inter- 
action strength U = 5.0 eV and J = 0.7eV^, and tem- 
perature T = 72.5 K. To explore the sensitivity of the 
magnetic moment to the strength of the Coulomb inter- 
action, we performed calculations for other values of U 
and J around the ab initio values. We found that the 
size of the magnetic moment can be well parameterized 
by the simple formula m = (0.4 U-\-7.2 J — 6.1ey)/iB/eV. 
Hence, magnetization is most sensitive to the value of the 
Hund's coupling J, rather than U. 
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The two types of polarization, spin and orbital, are 
dominant at different energy scales. Spin polarization 
is stronger at higher energy (far below the Fermi level) 
while orbital polarization is strong at low energy (close 
to the Fermi level). 

To demonstrate that the spin polarization is stronger 
at high energy, we plot in Fig. |3] the exchange splitting 
of iron 3d orbit als, defined by 

T^i,^{uj) - Si,;(a;). 

Here S(co') is the frequency dependent potential (self- 
energy) , which is added to the one particle Hamiltonian, 
to produce the many body spectra. The high frequency 
limit of the exchange splitting is on average three times 
of its zero frequency value. Hence, the magnetic moment 
comes mainly from the high energy region, where the ex- 
change splitting is large. 
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FIG. 3: Spin polarization of the frequency dependent self en- 
ergy (real part) Si,^(cj) — Si,|(u;) of the Fe 3dxz, 3dyz, Sdxy 
orbit als, and the average over all five 3d orbit als in the SDW 
state of BaFe2As2 calculated by LDA+DMFT. 



At energies within the SDW gap, the kinetic energy is 
dominant. To minimize the kinetic energy in the magnet- 
ically ordered state, electrons create a highway in the di- 
rection in which spins are antiferromagnetically ordered, 
while they remain slow in the ferromagnetic direction due 
to Pauli blocking. This gives rise to strong orbital po- 
larization of the iron 3d orbitals at energies close to the 
Fermi level. 

The orbital polarization is related to the frequency de- 
pendent Weiss field A(a;), which describes the hybridiza- 
tion of the iron atom with the rest of the system. In Fig. [4] 
we plot the difference of the hybridization between the 
xz and yz orbital, i.e., Axz{^) — Ayz{uj). In the para- 
magnetic state, the xz and yz hybridization are equal 
due to the tetragonal crystal structure, but they become 
different in the SDW state. As shown in Fig. [4] the or- 
bital polarization acquires finite value only in the region 
of dominant ly coherent spectra (within 1.5 eV from the 
Fermi level), and gets strongly enhanced within the SDW 
gap, where the quasiparticles are well defined. 
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FIG. 4: Orbital polarization of the frequency dependent hy- 
bridization Axz{(^) — Ayz{(jo) of the Fe 3dxz and 3dyz or- 
bitals summing over both spin channels in the SDW state 
of BaFe2As2 calculated by LDA+DMFT. 



